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Abstract 

Data recorded by the JADE experiment at the PETRA e^e^ cohider were used 
to measure the event shape observables thrust, heavy jet mass, wide and total jet 
broadening and the differential 2-jet rate in the Durham scheme. For the latter three 
observables, no experimental results have previously been presented at these energies. 
The distributions were compared with resummed QCD calulations (C'(ag)+NLLA), 
and the strong coupling constant as{Q) was determined at different energy scales 
Q = ^/s. The results, 

a, (22 GeV) = 0.16ll°;°}? , a, (35 GeV) = 0.143^°;°^^ , a, (44 GeV) = 0.1371°;°^? , 

are in agreement with previous combined results of PETRA albeit with smaller un- 
certainties. Together with corresponding data from LEP, the energy dependence of 
as is significantly tested and is found to be in good agreement with the QCD expec- 
tation. Similarly, mean values of the observables were compared to analytic QCD 
predictions where hadronisation effects are absorbed in calculable power corrections. 
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1 Introduction 



Summaries of measurements of as from various processes and at different energy scales Q 
demonstrate that the energy dependence of as{Q) is in good agreement with the pre- 
diction of Quantum Chromodynamics (QCD). The uncertainties of these measurements, 
both experimental and theoretical, are different and their correlations are, in general, not 
known P,^. More quantitative studies of the running of as therefore require the existence 
of consistent measurements over large ranges of the energy scale Q, for the same process, 
using identical experimental techniques and theoretical calculations in order to minimise 
point-to-point systematic uncertainties. 

Significant progress has been made in perturbative QCD calculations since 1992. Ob- 
servables have been proposed for which perturbative predictions are extended beyond the 
next-to-leading-order (0(q;^)) through the inclusion of leading and next-to-leading 
logarithms which are summed to all orders of as (NLLA) 1^-0 . These calculations exhibit 
a better stability to contributions of unknown higher order corrections, which are usually 
estimated by variations of the renormalisation scale /x. 

The experiments at LEP and SLC provided a number of significant determinations of 
as from hadronic event shapes and jet production, based on 0{a'^)+NLLA calculations, 
at centre-of-mass energies y/s at and above 91 GeV, the mass of the Z° boson. Detailed 
studies of the high statistics data samples from the LEP experiments provide a better 
understanding of the phenomenology of the hadronisation process and thus of the mod- 
elling of hadronic final states with Monte Carlo programs. Determinations of as at LEP 
and SLC []8Hl^ therefore have smaller uncertainties than those which are available from 
previous measurements at lower e"'"e~ centre-of-mass energies |]1],|1^. There are only a 
few recent measurements of as at lower energies. These either employed only some of the 
observables which are now available |T3,|TH| or they were based on limited samples of Z° 



decays with final state photon radiation |16|. Therefore equivalent studies with data in 



the centre-of-mass energy range from ^/s = 14 to 46.7 GeV, taken at the PETRA collider 
which was shut down in 1986, are desirable. 

In this paper we present an (9(a^)-|-NLLA determination of as at ^/s = 22, 35, and 
44 GeV using data from the JADE experiment [0,113 at PETRA. The selection of the 



JADE data and Monte Carlo event samples are described in Section The measurement 
of event shape distributions, the corrections for detector imperfections and for initial state 
photon radiation as well as the estimate of the experimental uncertainties are outlined in 
Section ^ The corrected event shape distributions and the determination of the strong 
coupling constant as{Q) are presented in Section |[ A study of the energy dependence 
of mean values of event shape distributions and their comparison with analytic QCD 
calculations comprising power corrections to account for hadronisation effects is presented 
in Section |. In Section ^ the results are summarised and the conclusions are drawn. 



2 Data samples and Monte Carlo simulation 

For the studies presented in this paper, we analysed data recorded with the JADE detector 
in 1981, 1984 to 1985, and 1986 at centre-of-mass energies of 22 GeV, 39.5-46.7 GeV, and 
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around 35 GeV, respectively. The JADE detector was one of the five experiments at the 
PETRA electron-positron collider. It was operated from 1979 until 1986 at centre-of-mass 
energies of y/s = 12 to 46.7 GeV. A detailed description of the JADE detector can be found 
in |]l],0. The main components of the detector were the central jet chamber to measure 
charged particle tracks and the lead glass calorimeter to measure energy depositions of 
electromagnetic showers, which both covered almost the whole solid angle of Att. 



Multihadronic events were selected by the standard JADE selection cuts [0 which were 
based on minimum energy deposits in the calorimeter and a minimum number of tracks 
emanating from the interaction region. All charged particle tracks with a total momentum 
of IpI > 100 MeV/c were considered in the analysis. Energy clusters in the electromagnetic 
calorimeter were considered if their energies exceeded 150 MeV after correction for energy 
deposited by associated tracks. Charged particle tracks were assumed to be pions while 
the photon hypothesis was assigned to electromagnetic energy clusters. 

In order to remove background from two-photon processes and r-pair events and from 
events which lost a substantial part of their energy due to hard initial state photon ra- 
diation, further constraints were imposed on the visible energy i?vis = the total 
missing momentum Pmiss = I Pi I {Pi and Ei are the 3-momentum and the energy of the 
tracks and clusters), the longitudinal balance relative to the e"'"e~ beam axis of momenta 
Pbai = I J2Pi/E^is\ and the polar angle of the thrust axis, 9t- 

• Evis > y/s/2 ; 

• Pmiss < 0.3 ■ ; 

• Pbai < 0.4 ; 

• I cos6'r| < 0.8 . 

With these cuts, the backgrounds from 77 and r-pair events were reduced to less than 
0.1% and 1%, respectively [|19|. The final numbers of events which were retained for this 
analysis are listed in Table 

The retrieval of data files eleven years after shutdown of the experiment was difficult 
and turned out to be incomplete at this stage of the analysis. Comparisons of the numbers 
given in Table |1] with previous JADE publications |T9|-[2T[] revealed that we were missing 
data sets of about 250 events around 22 GeV and about 450 events around 44 GeV. In 
addition, the original files containing information about the luminosity of different running 
periods could not be retrieved, so that only approximate values of integrated luminosities 



year 


[GeV] 


data 


MC 


1981 


22 


1404 




1984/85 


40-48 


6158 


14497 


1986 


35 


20 926 


25123 



Table 1: Number of events in data and in Monte Carlo detector simulation retained after 
application of the multihadron selection cuts described in the text. 
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Ref. [Tg] 


this analysis 


22 GeV 






R2 


72.5 ± 1.2 


72.7 ± 1.2 


R3 


27.1 ± 1.2 


27.0 ± 1.2 




0.42 ±0.16 


0.28 ±0.14 


35 GeV 






R2 


77.7 ±0.4 


78.2 ±0.3 


R3 


22.0 ±0.4 


21.6 ±0.3 


R4 


0.31 ±0.05 


0.24 ±0.03 


44 GeV 






R2 


79.8 ±0.5 


79.2 ±0.5 


R3 


20.1 ±0.5 


20.7 ±0.5 


i?4 


0.14 ±0.05 


0.15 ±0.05 



Table 2: A comparison of relative n-jet production rates as percentages of all hadronic 
events, using the JADE jet finding algorithm with y^ut = 0.08 
applied; the errors are statistical only. 



19 . No corrections are 



corresponding to our final number of events can be givenQ: the data samples shown in 
Tab. |I] correspond to about 2.4 ph-\ 80 pb'^ and 40 pb'^ at 22 GeV, 35 GeV and 
44 GeV centre-of-mass energy, respectively. 

In order to verify the compatibility of this study with results which were previously 
published by JADE, we repeated a determination of the relative 2-, 3- and 4-jet event 
production rates as published in |]19[, using the original JADE jet finder with a resolution 



parameter of i/cut = 0.08. The results are presented in Table |^. Considering the fact that 
our present data samples at 22 GeV and 44 GeV lack about 10% of the original ones and 
that the samples around 35 GeV are from different running periods (1986 for this analysis, 
1984-1985 for Reference [0), the agreement between the old and this new study is very 
good. This demonstrates that we are able to perform detailed studies of event properties 
in a consistent way. 

Corresponding Monte Carlo detector simulation data were retrieved for 35 and 44 GeV. 



They were generated using the QCD parton shower event generator JETSET 6.3 p2[. The 
Monte Carlo events at 35 GeV were generated using the coherent branching for the parton 
shower while the 44 GeV events had non-coherent branching |. The main parameters used 



for event generation are given in Section 3.3. Both samples included a simulation of the 



acceptance and resolution of the JADE detector. 

Comparisons of the measured and simulated distributions of visible energy, momentum 
balance, missing momentum and other quantities showed that the Monte Carlo simulation 
gave a reasonable description of the measurements. The simulated data can thus be 



^It is not expected that the missing events alter the measured distributions in a systematic manner. 
Also, detailed knowledge of luminosities is not required for the following studies of normalised event shape 
distributions. 

^The different treatment of coherence in these samples of simulated data has no visible influence on 
the results of this study; see also Fig. 1, Figs. 3-5 and the discussions in Section 0. 
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used to correct for detector effects in the measured data. As an example we show in 
Figure |l] the distributions of the thrust observable 1 — T and of the differential two-jet 
rates D21 measured at 35 and at 44 GeV. The definitions of these observables are given in 
Section |3.1| . In general, we found a good agreement of the detector simulation with data for 
all event shape distributions studied here, irrespective of coherent or non-coherent parton 
branching. 

3 Experimental procedure 

3.1 Event shapes and differential 2-jet rate 

From the data samples described in the previous section, the event shape distributions 
of thrust, the heavy jet mass, the total and wide jet broadening and the differential 2- 
jet event rate using the Durham jet finder were determined. For convenience we list the 
definitions of these observables. 



The vector n which maximises the expression in parentheses is the thrust axis n^. 
It is used to divide an event into two hemispheres Hi and H2 by a plane through 
the origin and perpendicular to the thrust axis. 

Heavy Jet Mass Mh'- 

From the particles in each of the two hemispheres defined by the thrust axis an 
invariant mass is calculated. The heavy jet mass Mh [0,^ is defined by the larger 
of the two masses. This analysis used the measured heavy jet mass scaled by the 
visible energy E^i^ which is, after correction for detector resolution, acceptance, and 
for initial state radiation, equal to Mh / \/s- 

Jet Broadening B: 

The jet broadening measures are calculated by the expression 0: 



for each of the two hemispheres, H^, defined above. The total jet broadening is given 
by = Bi + 82- The wide jet broadening is defined by = max(i?i, B2). 

Durham differential 2-jet rate D2'- 

Jets are reconstructed by a standard recombination algorithm: For any combination 
of two particles i and j in an event a measure of distance, yij, is calculated according 
to the Durham recombination scheme 



Thrust T: 



The thrust value of a hadronic event is defined by the expression |^ 





where Ei and Ej are the energies of the particles and cos 9ij is the angle between 
their 3- momentum vectors. The pair i, j oi particles with the smallest value of yij 
is replaced by a pseudoparticle k with 4-momentum = pi +Pj- This procedure 
is repeated until exactly three pseudoparticles remain which are called jets. The 
smallest i/ij corresponding to these three jets is indicated by 1/23 throughout the 
paper. At this particular value the number of reconstructed jets changes from 3 to 
2. D2 is the normalised differential cross-section as a function of ^23 [0- 

In the following we use the symbols T, M^, Bt, Bw and D2 to denote thrust, heavy jet 
mass, total and wide jet broadening, and the differential 2-jet rate, respectively. 

3.2 Correction procedure 

The event shape data were corrected for the limited acceptance and resolution of the 
detector and for initial state photon radiation effects by applying a bin-by-bin correction 
procedure. Correction factors were defined by the ratio of the distribution calculated from 
events generated by JETSET 6.3 at hadron level over the same distribution at detector 
level. The hadron level distributions were obtained from JETSET 6.3 generator runs 
without detector simulation and without initial state radiation, using all particles with 
lifetimes r > 3 ■ 10"^'' s. The model events at detector level contained initial state photon 
radiation and a detailed simulation of the detector response and were processed in the 
same way as the data. 

In a second step, the data distributions were further corrected for hadronisation ef- 
fects. This was done by applying bin-by-bin correction factors derived from the ratio of 
the distribution at parton level over the same distribution at hadron level, which were cal- 
culated from JETSET generated events before and after hadronisation, respectively. The 
data distributions, thus corrected to the parton level, can be compared to analytic QCD 
calculations. 



3.3 Systematic uncertainties 

To study systematic uncertainties of the corrected data distributions we modified details of 
the event selection and of the correction procedure. For each variation the whole analysis 
was repeated and any deviation from the main result was considered a systematic error. 
In general, the maximum deviation from the main result for each kind of variation was 
regarded as symmetric systematic uncertainty. The main result was obtained using the 
default selection and correction procedure as described above. 

We restricted the measurement of the event shape distributions to rely either on tracks 
or on clusters only. We varied the cut on cosOt by ±0.1. The cut on Pmiss was either 
removed or tightened to Pmiss < 0.25 ■ ^/s. Similarly, the momentum balance requirement 
was either restricted to pbai < 0.3 or dropped. We also varied the cut for the visible energy 
E^is by ±0.05 • ^/s. In order to check the residual contributions from r-pair events we also 
required at least seven well-measured charged tracks. 

To study the impact of the hadronisation model of the JETSET 6.3 generator, the 
values of several significant model parameters were varied around their tuned default 



values used for our main result. We took the tuned values from Reference 21 , where 



6 



the QCD parameter was Alla = 400 MeV, the cut-off for the parton shower development 
was Qq = 1 GeV, and the width of the transverse momentum distribution of the hadrons 
with respect to the direction of the quark was ctq = 300 MeV. According to this Reference 
we chose the LUND symmetric fragmentation function with a = 0.5 and b = 0.9 for the 
fragmentation of the hght u, d, and s quarks. The heavy c and b quark were fragmented 
applying the Peterson et al. |28| fragmentation function using ec = 0.05 and eb = 0.01 |21|. 

Different sets of correction factors to correct the data from hadron level to parton level 
were generated by varying single parameters of the JETSET generator. The variations 
were chosen to be similar to the one standard deviation percentage limits obtained by the 
OPAL Collaboration from a parameter tuning of JETSET at ^/s = M^o . 



In particular, we investigated the effects due to parton shower, hadronisation param- 
eters, and quark masses. The amount of gluon radiation during the parton shower devel- 
opment was modified by varying Alla by ±50 MeV. To vary the onset of hadronisation, 
we altered the parton shower cut-off parameter Qq by ±0.5 GeV. We used the full ob- 
served variation of to reflect a variation of Qq between and 2 GeV. The width of the 
transverse momentum distribution in the hadronisation process was varied by ±30 MeV. 
The LUND symmetric fragmentation function was varied by changing the a parameter by 
±0.225 whereas the h parameter was kept fixed. As a systematic variation we used the 
LUND fragmentation function also for charm and bottom quarks. The effects due to the 
bottom quark mass were studied by restricting the model calculations which were used to 
determine the correction factors to up, down, strange, and charm quarks (udsc) only. In 
this case, any deviation from our main result was treated as asymmetric error. 

No detector level Monte Carlo simulation data were available for the 22 GeV data. In 
order to obtain consistent detector corrections also at this energy, we studied the energy 
dependence of the detector correction from the 35 and 44 GeV Monte Carlo samples. Here 
we considered only the differential 2-jet rate, D21 of the Durham jet finder scheme because 
it is known to depend to a lesser extent on hadronisation and detector effects. In Figure |^ 
the detector correction factors as obtained at 35 GeV and at 44 GeV are displayed. In 
general, the corrections are small, and their size is about the same at both centre-of-mass 
energies. There is no apparent energy dependence of the detector correction within the 
range, indicated by the arrow, which was considered for the fit of as at 22 GeV. We 
therefore applied the 35 GeV detector correction to the differential 2-jet rate measured 
from the 22 GeV data, and studied only the well known dominating sources of systematic 
uncertainties. The correction of hadronisation effects was then determined from JETSET 
generator runs at 22 GeV centre-of-mass energies, exactly as for the data at 35 and 44 GeV. 



4 Determination of ag 

4.1 Corrected event shape distributions 

After applying the corrections for detector and for initial state radiation effects we ob- 
tained the event shape distributions at hadron level. These are shown in Figures ^, ^ and 
^ for the 44, 35, and 22 GeV data samples. For comparison the respective distributions 
predicted by the JETSET 6.3 generator, at hadron level, are also shown. There is excel- 
lent agreement between the data and the model over the whole kinematic range of the 
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observables. In Tables ^, | and ^ the corrected data values are listed with statistical errors 
and experimental systematic uncertainties. The mean values of the distributions are also 
given. 



4.2 QCD calculations for event shapes 

The distributions of the event shape observables used in this analysis are predicted in 
perturbative QCD by a combination of the Ola"^) 0] and the NLLA ^-0] calculations. 
The 0{al) calculation yields an expression of the form 



Roiady) = ^ + My) {^)+Biy) (^)' 



where R{y) = Jq dy l/aQ-da/dy is the cumulative cross-section of an event shape observ- 
able y normalised to the lowest order Born cross-section ctq. The NLLA calculations give 
an expression for R{y) in the form: 



a 



RNLLAiy) = U + ^1 ( ^ j + ^4^) ) • 



2^ 



where L = ln(l/y). The functions gi and g2 are given by the NLLA calculations. The 
coefficients Ci and C2 are known from the 0{a'^) matrix elements. 



4.3 Determination of ag using 0(ag)+NLLA calculations 

We determined by fits to event shape distributions of 1 — T, Mh, Bt, and of D2 
corrected to the parton level. For the sake of direct comparison to other published results 



we closely followed the procedures described in |[TT|,pO|,RT[]. We chose the so-called ln(i?)- 



matching scheme to merge the 0{al) with the NLLA calculations. The renormalisation 
scale factor, x^ = ^/ y/s, was set to = 1 for what we chose to be the main result. Here, 
the value of /i defines the energy scale at which the theory is renormalised. 

The fit ranges for each observable were determined by choosing the largest range for 
which the hadronisation uncertainties remained below about 10 %, for which the x^/d.o.f. 
of the fits did not exceed the minimum by more than a factor of two, and by aiming 
at results for that are independent of the fit range. The remaining changes when 
enlarging or reducing the fit range by one bin on either side were taken as systematic 
uncertainties. Only statistical errors were considered in the fit thus resulting in x^/d.o.f. 
larger than unity. The finally selected fit ranges, the results of the fits and of the study 
of systematic uncertainties are tabulated in Tables |^, ||, and |^ and are shown in Figures ^, 
0, and g 

The dependence of the fit result for as on indicates the importance of higher order 
terms in the theory. We also changed the renormalisation scale factor in the range of 
Xfj^ = 0.5 to 2.0. We found variations of similar size as the uncertainties from the detector 
correction and the hadronisation model dependence. The differential 2-jet rate, D2, in the 
Durham jet scheme exhibits the smallest renormalisation scale uncertainties, resulting in 
the smallest total error of all observables considered in this analysis. The values of as and 
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the errors obtained at 35 and 44 GeV are shown in Figure 0. In these diagrams also the 



values measured by the OPAL Collaboration at ^/s = Mgo [|ri[] are shown for comparison. 
The values of exhibit a similar scattering pattern at all energies. This demonstrates 
the strong correlation of the systematic uncertainties, which are dominated by theoretical 
and hadronisation uncertainties. 

The individual results of the four event shape observables and the differential 2-jet rate 



were combined into a single value following the procedure described in References 



32| . This procedure accounts for correlations of the systematic uncertainties. At each 



energy, a weighted average of the five values was calculated with the reciprocal of the 
square of the respective total error used as a weight. In the case of asymmetric errors 
we took the average of the positive and negative error to determine the weight. For each 
of the systematic checks, the mean of the ag values from all considered observables was 
determined. Any deviation of this mean from the weighted average of the main result was 
taken as a systematic uncertainty. 

With this procedure we obtained as final results for ag 

as{UGeV) = 0.1372 ± 0.0017(stat.) l°;°J°^(syst.) 
a, (35 GeV) = 0.1434 ± 0.0010(stat.) l[j:[|J^^(syst.) 
a, (22 GeV) = 0.1608 ± 0.0083(stat.) +°;°J^9(gyg^_) ^ 

where the result at 22 GeV is based on the differential 2-jet rate only. The systematic errors 
at 44, 35, and 22 GeV are the quadratic sums of the experimental uncertainties (±0.0034 
±0.0018, ±0.0030), the effects due to the Monte Carlo modelling (tHo2?, -omti -omit 
and the contributions due to the variation of the renormalisation scale I'+^-o*'^^ +o.oo87 



0.0027) -0.0034) -0.0056) 
0.0082 +0.0087 
0.0054) -0.0052) 

"aoool)- It should be noted that the modelling uncertainties due to quark mass effects 



contribute significantly to the total error. 



4.4 Determination of ag using 0{a'^) calculations 

For comparison, we repeated the fits using fixed order 0{al) calculations only. The fit 
ranges for each distribution had to be readjusted in order to match the stability require- 
ments given aboveQ. All systematic checks were done as described above except for the 
variation of the renormalisation scale factor x^. Instead, the (9(a^)-fits were performed 
once with fixed to 1 and once with as a free parameter of the fit. The fit ranges 
were the same in both fits except for D2 where it had to be enlarged towards the lower 
end in order to obtain a stable fit with free parameter. The mean value of as from 

the two fits was taken as the final result while half of the difference between the two was 
assigned as a systematic error due to the unknown higher orders in perturbation theory. 
The results of the 0{a1) fits are summarised in Tables ^ and |TT]. The corresponding 
theoretical predictions were superimposed on the results of the (9(af)+NLLA fits that are 
presented in Figures ^ and |^. All results at a given centre-of-mass energy agree with each 

•^From corresponding studies at LEP ||ll|, ^ it is known that different fit ranges are required for the 
0{ag) and for C'(a^)+NLLA predictions. This is also supported by theoretical considerations, since the 
inclusion of NLLA is supposed to extend the degree of reliability especially in the 2-jet region of phase 
space, i.e. at small values of the event shape observables used in this study. 
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other but the values from the (9(a^)+NLLA fits are systematically lower. Again, the 
pattern between these results and those obtained at the higher energies of LEP |^ is very 
similar. 



5 Mean Values of Distributions and QCD Power Cor- 
rections 

5.1 Power corrections 

The value of as can also be assessed by the energy dependence of mean values of event 
shape distributions. Presently, the mean values of the observables considered in this 
analysis are calculated up to 0{a1). For an observable the perturbative prediction is 



where the coefficients Aj: and Bj: were determined from the OipP'^ perturbative calcula- 
tions [m^,|2^,|3^. The term —2Aj: accounts for the difference between the total cross- 
section used in the measurement and the Born level cross-section used in the perturbative 
calculation. The numerical values of these coefficients are summarised in Table ^ 

Instead of correcting for hadronisation effects with a Monte Carlo event generator as 
we did for the (Xg determination presented in Section ^ we considered additive power- 
suppressed corrections (1/ (a/s)^) to the perturbative predictions of the mean values of the 
event shape observables. Such corrections are expected on general grounds for hadro- 
nisation and other non-perturbative effects, for example renormalons [0. The non- 



perturbative effects are due to the emission of very low energetic gluons which can not 
be treated perturbatively due to the divergence of the perturbative expressions for at 
low scales. In the calculations of Reference which we used in this analysis a non- 
perturbative parameter 

- ( \ 

was introduced to replace the divergent portion of the perturbative expression for ^^(y^) 
below an infrared matching scale fij. The general form of the power correction to the 



dk as{k)-kP 



Observable JF 



A 



B 



P 



(T) 
(Mjj/s) 
(Bt) 
{Bw) 
(1/23) 



2.103 
2.103 
4.066 
4.066 
0.895 



44.99 
23.24 
64.24 
-9.53 
12.68 



-1 
1.0 ±0.5 
1.0 ±0.5 
1.0 ±0.5 





1 
1 



Table 3: Coefficients of the perturbative prediction [Q]^, Q and coefficients and 
parameters of the power corrections ||3^ to the mean values of the event shape observables. 
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mean value of an observable T assumes the form 

where Cp = 4/3. The factor Pq = (IIC^ — 2Nf )/3 stems from the QCD /5-function 
of the renormalisation group equation. It depends on the number of colours, Ca = 3, 
and number of active quark flavours Nf, for which we used Nj = 5 throughout the 
analysis. The term K = (67/18 — tt^/Q)Ca — 5/9 ■ A'^^ originates from the choice of 
the MS renormalisation scheme. The remaining coefficient ajr and the parameters p and r 
depend on the event shape observable. For completeness, these coefficients and parameters 
obtained in Reference |^ are also listed in Table 

5.2 Determination of as using power corrections 

We determined ^^(Mzo) by fits of the expression 

to the mean values of the five observables investigated in this analysis^ including the 
measured mean values obtained by other experiments at different centre-of-mass ener- 
gies 0,|3l|,^,^. For the central values of as from the fits we chose a renormalisation 



scale factor of a;^ = 1 and an infrared scale of yU/ = 2 GeV. The x^/d.o.f. of all fits were 
between 0.8 ((M^/s)) and 4.2 {{Bt)). We estimated the systematic uncertainties by vary- 
ing from 0.5 to 2 and from 1 to 3 GeV. Since the precision of the coefficient ajr as 



given in Reference |35] for the heavy jet mass Mjj and the two jet broadening measures, 
Bt and Bw is only ±50%, we assigned an additional uncertainty to as due to the variation 
of these coefficients by this amount. 

In the case of (?/23) no coefficient ajr is given in Reference We investigated the size 
of ajr by fitting with fixed to the world average aJ-^-{Mzo) = 0.118. All fits with 
p = 1 or 2, and r = or 1 resulted in very small values of ajr compatible with zero. From 
this we conclude that power corrections to the perturbative prediction for (1/23) can be 
neglected for the energy range considered. Therefore, we used the perturbative prediction 
only, for which we obtained a good fit with x^/d.o.f. = 1. 

The results of the fits are shown in Figure |10| and the numeric values are listed in 



Table |I2|. It presents the values for as and for uq, the experimental errors and systematic 
uncertainties of the fit results. We consider these results based on power corrections as 
a test of this theoretical prediction. It should be noted that the theoretically expected 



universality of ao is not observed. The issue of universality is further addressed in . 

Employing the procedure used in Section | to combine the individual as values, we 
obtained 

a,(Mzo) =0.1155 t°;0°^2 



^Our results for (M'fj/s) are 0.0745 ± 0.0011 and 0.0679 ± 0.0008 at = 35 GeV and 44 GeV 
respectively. The errors are the statistical and systematical uncertainties added in quadrature. 
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where the error is the experimental uncertainty (±0.0013), the renormahsation scale un- 
certainty (laooal), the uncertainty due to the choice of the infrared scale (laoo^g) and 
the uncertainties of the non-perturbative coefficients ajr (lo;oo2o)) combined in quadra- 
ture. This result is in good agreement with the world average value of a^'^^Mio) = 
0.118 ±0.006. Our value is also in agreement with the results of similar studies for different 



sets of observables by the DELPHI Collaboration |37] and by the HI Collaboration |39 



6 Summary and Conclusions 

Data recorded by the JADE experiment at centre-of-mass energies around 22, 35, and 
44 GeV were analysed in terms of event shape distributions and differential 2-jet rates. For 
most of the observables no experimental results have previously been presented, because 
the total and wide jet broadening, Bj- and B\y, as well as the Durham jet finding scheme 
were proposed only after the shutdown of the experiments at the PETRA accelerator. 

The measured distributions were corrected for detector and initial state photon radia- 
tion effects using original Monte Carlo simulation data for 35 and 44 GeV. The simulated 
data are based on the JETSET parton shower generator version 6.3. The same event gen- 
erator was also employed to correct the data for hadronisation effects in order to determine 
the strong coupling constant as- 

Our measurements of as are based on the most complete theoretical calculations avail- 
able to date. For all observables theoretical calculations exist in 0{al) and in the next- 
to-leading log approximation. These two calculations were combined using the ln(_R)- 
matching scheme. 

The final values of ag at the three different centre-of-mass energies are 

«,(44GeV) = 0.137^^0? 
a, (35 GeV) = 0.143 tofol 



a,(22 GeV) = 0.161 



-0.016 
-0.011 



where the errors are statistical, experimental systematics, Monte Carlo modelling and 
higher order QCD uncertainties added in quadrature. The dominant contributions to the 
total error came from the choice of the renormahsation scale and from uncertainties due 
to quark mass effects. 

The as result at 22 GeV was obtained from the differential 2-jet rate only. Note, 
however, that for 35 and 44 GeV the as value obtained from the differential 2-jet rate has 
the smallest total error and is very close to the weighted average as can be inferred from 
Figure |^. We therefore consider the value obtained at 22 GeV a good approximation 
of the projected result of a more comprehensive study at this energy. 

The fits for as were also performed using the Ola"^) calculation alone. All results were 
found to be consistent with each other. 

These results agree well with those which are available from previous measurements 
of as in the PETRA and PEP energy range; see e.g. [Q,[T3| for reviews of that time. Our 



results, however, include more detailed systematic studies, are based on more observables 
and use more advanced theoretical calculations; nevertheless they exhibit smaller total 
errors. 
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Similarities between the main components of the JADE detector at PETRA and 



the OPAL detector |Q at LEP, as well as between this analysis and studies performed by 



the OPAL Collaboration |]TT],^,|31]] at ^/s = 91.2, 133, and 161 GeV, suggest the energy 
dependence of as in the centre-of-mass energy range of ^/s = 22-161 GeV can be reliably 
tested, because the systematic uncertainties of these measurements are partly correlated. 

The as results from OPAL and from this analysis are shown in Figure |lT[ The result 
of a fit of the 0{al) QCD prediction to the data is shown by the solid line. The 
fit resulted in a,(Mzo) = 0.1207 ± 0.0012 and xV^-cf. = 4.9/5, taking into account 



only statistical and experimental uncertainties, which are displayed in Fig. |rT| as the solid, 
innermost error bars. The other systematic uncertainties, due to hadronisation and to 
unknown higher order contributions, are assumed to be fully correlated at all energies and 
thus are not considered in this test of the energy dependence of Og. A visible trend of the 
lower energy results all lying above and the higher energy ones lying below the fitted QCD 
curve can be consistently explained within the assigned experimental uncertainties which 
is indicated by the value of x^/d.o.f. = 1. 

A fit for the hypothesis of a constant value of as gives as = 0.1328 ± 0.0014 and 
X^/d.o.f. = 101/5, which has a vanishing probability. The energy dependence of as is 
therefore significantly demonstrated by the results from the combined JADE and OPAL 
data. 

Evolving our as measurements to y/s = Mzo the results obtained at 44, 35 and 22 GeV 
transform to 0.122 Iggg^, 0.122 and 0.124 respectively. The combination of 

these values gives ^^(Mzo) = 0.122 ^o;oo5. This value is consistent with the direct mea- 
surement at ^/s = Mzo by the OPAL Collaboration of as{Mzo) = 0.117t{][J{j^ for the 
same subset of observables. 

The energy dependence of the mean values of the distributions can be directly compared 



with analytic QCD predictions plus power corrections for hadronisation effects I^Sj. Until 
recently, such studies were hardly possible since for most of the observables no results were 
available at energies below the Z" mass scale. With the inclusion of the results presented 
in this paper, comprehensive fits of the analytic predictions to the data are now possible. 
Our studies resulted in 

a.(Mzo)= 0.116 1^05 

which is in good agreement with our results from the (9(a^)-|-NLLA fits, with measure- 



ments at LEP and at HERA ||39| and also with the world average value. 

In summary, new studies of hadronic final states of e"'"e~-annihilations in the PETRA 
energy range provided valuable information which was not available before. New results 
of as, obtained in a similar manner as those from the experiments at LEP, provide a 
significant test of the running of and thus of the non-abelian nature of QCD. Evolved 
to the Z° mass scale, the results are in good agreement with those obtained at LEP, and are 
of similar precision. A direct comparison of the energy dependence of the mean values of 
the measured distributions with analytic QCD calculations plus power corrections provide 
alternative ways to test QCD, without the need to rely on phenomenological hadronisation 
models. 



■''This value of as(Mzo) corresponds to a QCD scale ALI = 242 ± 15 MeV for five active quark flavours. 
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Work has been started to further decrease the overall uncertainties of the results pre- 
sented in this paper, and to study more aspects of QCD using the JADE data samples. 
This can be achieved by the use of more recent event generators and the JADE detec- 
tor simulation software. This will provide the possibility to study the data at the lowest 
PETRA energies, around ^/s = 14 and 22 GeV, in more detail, i.e. for energy scales at 
which the variation of is strongest. In addition, the significance of results from data 
at PETRA energies will increase from a better and more fundamental treatment of the 
b-quark mass, in theory E3] as well as in experiment. 
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Table 4: Event shape data at ^/s = 44 GeV for the observablcs described in the text. The 
values were corrected for detector and for initial state radiation effects. The first errors 
denote the statistical and the second the experimental systematic uncertainties. 



18 



1 

1 


-T 


l/a. 


da/d(l 




J- ) 




0.00- 


-0 


02 


0.638 


± 





028 


± 


0.181 


Mh 




1/cr 


■ da/d(MH/v^) 


0.02- 


-0 


04 


6.43 


± 





12 


± 


0.23 


0.00- 


-0.06 


0.002 


± 





000 


± 


0.004 


0.04- 


-0 


06 


11.00 


± 





17 


1 


0.25 


0.06- 


-0.10 


0.017 


± 





001 


± 


0.009 


O.Od- 


-0 


08 


9.47 


± 





16 


1 


0.43 


0.10- 


-0.14 


0.288 


± 





009 


± 


0.041 


0.08- 


-0 


10 


6.43 


± 





13 


1 


0.21 


0.14- 


-0.18 


2.566 


± 





043 


± 


0.095 


0.10- 


-0 


12 


4.049 


± 





095 


1 


0.173 


0.18- 


-0.22 


6.278 


± 





090 


± 


0.361 


0.12- 


-0 


14 


3.033 


± 





084 


1 


0.239 


0.22- 


-0.26 


5.463 


± 





088 


± 


0.319 


0.14- 


-0 


16 


1.962 


± 





065 


1 


0.151 


0.26- 


-0.30 


3.823 


± 





073 


± 


0.173 


0.16- 


-0 


18 


1.704 


± 





062 


1 


0.179 


0.30- 


-0.34 


2.390 


± 





056 


± 


0.084 


r\ "1 o 

0.18- 


-0 


20 


1.209 


± 





050 


1 


0.117 


0.34- 


-0.38 


1.643 


1 

± 


u 


U4/ 


1 

± 


U.U88 


0.20- 


-0 


23 


n ooo 


1 


n 
u 


Uoo 


1 

± 


0.066 


U.oo- 


n /1 9 
-U.4Z 


1 nns 
i.UUo 


± 





037 


± 


0.039 


n 9^^ 


-0 


26 


U. ( 04 


_i_ 
It 


n 
u 


(Tie; 
Uoo 


-1- 


u.uyo 


n /19 


n AP. 

-U.4D 


0.626 


± 





030 


± 


0.040 


96- 


-0 


30 


453 


± 





022 


1 


049 


0.46- 


-0.50 


0.427 


lb 





025 


± 


0.067 


0.30- 


-0 


35 


0.186 


± 





012 


± 


0.061 


0.50- 


-0.55 


0.153 


± 





013 


± 


0.034 


0.35- 


-0.40 


0.070 


± 





008 


± 


0.010 


0.55- 


-0.60 


0.020 


± 





004 


± 


0.005 


0.40- 


-0.50 


0.001 


± 





001 


± 


0.001 


mean 


value 


0.2601 ± 0.0006 ± 0.0016 


mean 


value 


0.0938 


± 0.0004 ± 0.0015 





B 


T 




l/a- 








0.00- 


-0 


03 


0.018 


± 





003 


± 





028 


B 


w 




l/a. 


da/dBw 


0.03- 


-0 


06 


0.381 


± 





019 


± 





103 


0.00- 


-0 


02 


0.041 


± 





006 


± 


0.046 


0.06- 


-0 


08 


3.371 


± 





089 


± 





279 


0.02- 


-0 


04 


1.628 


± 





059 


± 


0.276 


0.08- 


-0 


10 


8.02 


± 





15 


± 





82 


0.04- 


-0 


06 


11.79 


± 





18 


± 


0.46 


0.10- 


-0 


12 


8.50 


± 





15 


± 





16 


0.06- 


-0 


08 


12.18 


± 





17 


± 


0.31 


0.12- 


-0 


14 


7.38 


± 





14 


± 





31 


0.08- 


-0 


10 


8.87 


± 





14 


± 


0.42 


0.14- 


-0 


16 


6.27 


± 





12 


± 





33 


0.10- 


-0 


12 


5.11 


± 





10 


± 


0.20 


0.16- 


-0 


18 


4.52 


± 





10 


± 





13 


0.12- 


-0 


14 


3.63 


± 





088 


± 


0.255 


0.18- 


-0 


20 


3.267 


± 





084 


± 





149 


0.14- 


-0 


16 


2.479 


± 





074 


± 


0.184 


0.20- 


-0 


22 


2.429 


± 





072 


± 





202 


0.16- 


-0 


18 


1.631 


± 





059 


± 


0.286 


0.22- 


-0 


24 


1.748 


± 





060 


± 





149 


0.18- 


-0 


20 


1.092 


± 





049 


± 


0.052 


0.24- 


-0 


27 


1.277 


± 





042 


± 





090 


0.20- 


-0 


23 


0.739 


± 





035 


± 


0.133 


0.27- 


-0 


30 


0.811 


± 





033 


± 





061 


0.23- 


-0 


26 


0.276 


± 





021 


± 


0.051 


0.30- 


■0 


35 


0.262 


± 





013 


± 





047 


0.26- 


-0 


30 


0.020 


± 





004 


± 


0.004 


0.35- 


-0.40 


0.020 


± 





003 


± 





006 


mean 


value 


0.0906 ± 0.0003 ± 0.0009 


mean 


value 


0.1439 ± 0.0004 ± 0.0012 





Table 5: Event shape data as for Table ^ but measured at ^/s = 35 GeV. 
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1/a ■ da/dy23 


OO KjrC V 


1/a ■ da/dy23 


0.000-0.001 


15.2 ± 


1.1 


± 


3.7 


0.000-0.001 


8.45 ± 


0.47 


± 


3 


23 


0.001-0.002 


67.4 ± 


3.1 


± 


9.0 


A AA1 A AAA 

0.001-0.002 


46.5 ± 


1.4 


± 


3 


6 


0.002-0.005 


86.7 ± 


2.3 


± 


4.1 


A AAA A AAr 

0.002-0.005 


73.5 ± 


1.1 


± 


2 


7 


A A A CT A A 1 A 

U. 005-0. 010 


49.4 ± 


1.4 


± 


5.9 


A AAC A A1 A 

0.005-0.010 


45.74 ± 


0.67 


± 


2 


13 


A A1 A A AAA 

0.010-0.020 


16.00 ± 


0.51 


± 


0.98 


A A 1 A A AAA 

0.010-0.020 


19.32 ± 


0.30 


± 





66 


A AAA A AOA 

0.020-0.030 


7.26 ± 


0.32 


± 


0.97 


A AAA A AOA 

0.020-0.030 


8.09 ± 


0.19 


± 





16 


A AOA A A/1A 

0.030-0.040 


4.04 ± 


0.23 


± 


0.69 


A AOA AA/1A 

0.030-0.040 


4.75 ± 


0.15 


± 





30 


A A/1A A AP'A 

0.040-0.050 


2.35 ± 


0.18 


± 


0.48 


AA/1A AAP'A 

0.040-0.050 


3.09 ± 


0.12 


± 





19 


A Ar'A A A/^A 

0.050-0.060 


1.94 ± 


0.16 


± 


0.61 


A A A A A/^A 

0.050-0.060 


2.39 ± 


0.11 


± 





14 


A A/^A A AOA 

0.060-0.080 


1.302 ± 


0.094 


± 


0.078 


A A/^A A AOA 

0.060-0.080 


1.568 ± 


0.060 


± 





098 


A AOA A 1 AA 

0.080-0.100 


0.897 ± 


0.079 


± 


0.081 


A AOA A 1 AA 

0.080-0.100 


1.096 ± 


0.052 


± 





109 


A 1 AA A 1 OA 

0.100-0.130 


U.6U1 ± 


U.Uo3 


1 

± 


n n f; /I 
U.Uo4 


A 1 AA A 1 OA 

0.100-0.130 


U.ol6 ± 


U.Uoo 


1 

± 


u 


1/11 


n 1 '^n n 1 fin 


0.484 ± 


0.053 


± 


0.171 


n T^n n 1 fin 


0.449 ± 


0.027 


± 





048 


160-0 200 


0.276 ± 


0.035 


± 


0.034 


160-0 200 


0.365 ± 


0.023 


± 





083 


0.200-0.250 


0.142 ± 


0.022 


± 


0.076 


0.200-0.250 


0.170 ± 


0.014 


± 





039 


0.250-0.400 


0.024 ± 


0.006 


± 


0.023 


0.250-0.400 


0.026 ± 


0.004 


± 





005 


mean value 


0.0229 ± 0.0005 ± 0.0019 


mean value 


0.0266 ± 0.0003 ± 0.0015 



22 GeV 


1/a- da/dy23 


0.000-0.001 


1.91 ± 


0.87 


± 


2.84 


0.001-0.002 


18.5 ± 


3.6 


± 


7.7 


0.002-0.005 


37.4 ± 


3.2 


± 


10.7 


0.005-0.010 


41.4 ± 


2.5 


± 


1.6 


0.010-0.020 


24.4 ± 


1.4 


± 


3.5 


0.020-0.030 


11.64 ± 


0.90 


± 


0.12 


0.030-0.040 


7.61 ± 


0.77 


± 


1.57 


0.040-0.050 


3.76 ± 


0.51 


± 


1.27 


0.050-0.060 


3.206 ± 


0.490 


± 


0.086 


0.060-0.080 


2.38 ± 


0.30 


± 


0.38 


0.080-0.100 


1.17 ± 


0.21 


± 


0.81 


0.100-0.130 


0.85 ± 


0.14 


± 


0.19 


0.130-0.160 


0.49 ± 


0.11 


± 


0.11 


0.160-0.200 


0.278 ± 


0.077 


± 


0.070 


0.200-0.250 


0.205 ± 


0.062 


± 


0.066 


0.250-0.400 


0.022 ± 


0.016 


± 


0.042 


mean value 


0.0311 ± 0.0011 ± 0.0018 



Table 6: Differential 2-jet rate D2 at ^ = 44 GeV, at 35 GeV and at 22 GeV. The 
values were corrected for detector and for initial state radiation effects. The first errors 
denote the statistical and the second the experimental systematic uncertainties. 
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1-T 


Mh 


Bt 


By/ 


D2 


a, (44 GeV) 


0.1457 


0.1423 


0.1417 


0.1278 


0.1344 


fit range 


0.08-0.3 


0.22-0.46 


0.080-0.27 


0.06-0.16 


0.005-0.200 


xVd-o.f. 


9.4/8 


31.0/5 


27.1/8 


48.8/4 


37.8/10 


Statistical error 


±0.0017 


±0.0017 


±0.0014 


±0.0016 


±0.0019 


lidL^ivo Uiiij 


U.UUl 1 


—0 001 9 


—0 0097 

U.UUZ, 1 


—0 0004 
u.uuuo 


—0 0099 
u.uuzz 




-i-D nn97 




-1-0 OOOQ 

n^U.UUUc/ 


-1-0 001 5 
n^u.uuxo 


-1-0 001 (S 
n^u.uu ±u 




+n nnni 

znu.uuux 


+n nnnp 


+n nnn"! 


+n nnns 

znu.uuuo 


+n nnnfi 

znu.uuuu 


Pmiss 


-i-n nnnfi 

znu.uuuu 


-I-n nnns 


-I-n nnni 


-I-n nnn9 

znu.uuuz 


-I-n nnn"^ 

znu.uuuu 


T)u 1 


±0 0004 

• yjyjyj^ 


±0 0004 

—1—\J •kJkJkJ^ 


±0 0003 


±0 0004 


±0 0002 

_l_ w • w w w ^ 


Nch 


+0.0002 


+0.0004 


+0.0003 


+0.0003 


+0.0006 


E ■ 

-'-'VIS 


±0.0003 


±0.0001 


±0.0003 


±0.0002 


±0.0002 


fit range 


±0.0015 


±0.0032 


±0.0023 


±0.0048 


±0.0032 


Experimental syst. 


±0.0032 


±0.0050 


±0.0037 


±0.0051 


±0.0040 




-1-0 nni Q 


-I-n nni r 


-I-n nni 7 


-I-n nni n 

n^u.uuxu 


n nnni 

\ u.uuux 




— U.UUIU 


— n nn9i 


— n nni 7 

— U.UUl / 


— n nnnQ 

— u.uuuy 


— n nnnii 

— u.uuuo 


rr — "^n A/TpV 


-1-0 nni n 


-I-n nnni 


-I-n nnnn 


-I-n nnn7 


-I-n nnn/1 

n^u.uuu^ 


rr 4- "^C] A/TpV 


— n nnnq 


— n nnn'? 

u.uuuo 


— n nni 1 

U .UU -L -L 


— n nnnfi 

u.uuuu 


— n nnno 

u.uuuy 


1j IN oV iiiilidi 


-I-n nm 9 


-I-n nnnq 

n^U.UUUc/ 


-i-n nni 7 

n^U.UU-L / 


-I-n nni 5 

n^U.UU-LtJ 


-I-n nnn^ 

n^u.uuuo 


icj/Q OUU IVit: V 


— n nnns 


-I-n nni i 


— n nnn7 

— U.UUU ( 


-I-n nni 9 

n^u.uuiz 


-I-n nn9fi 

n^U.UUZL) 




1(1 000 Q 

\\J .\J\J\JO 


0007 


1 0009 


0009 


001 5 


A _ 50 MpV 


—0 0005 


-1-0 0001 


—0 0013 


-1-0 0001 

1 \j . www X 


-1-0 0003 

\^ W . W W WO 


A + 50 MeV 


+0.0008 


-0.0005 


+0.0008 


< 0.0001 


-0.0011 


udsc only 


+0.0040 


+0.0007 


+0.0064 


+0.0047 


+0.0049 


ML statistics 


±0.0011 


±0.0011 


±0.0009 


±0.0011 


±0.0012 


MC modelling 


+U.UUoU 
-0.0030 


+U.UU00 
-0.0032 


1 n nn'70 
+U.UU/ z 

-0.0032 


+U.UUo4 
-0.0027 


+U.UUO ( 
-0.0045 


= 0.5 


—0.0089 


— 0.00o7 


—0.0100 


—0.0065 


—0.0007 


a;^ = 2.0 


+0.0115 


+0.0092 


+0.0125 


+0.0082 


+0.0045 


Total error 


+0.0131 
-0.0101 


+0.0111 
-0.0091 


+0.0150 
-0.0112 


+0.0112 
-0.0088 


+0.0091 
-0.0063 



Table 7: Values of ci;5(44 GeV) derived using the C(Q;f)+NLLA QCD calculations with 
X/^ — 1 and the ln(i?)-matching scheme, fit ranges and x^/d-O.f. values for each of the 

five event shape observables. In addition, the statistical and systematic uncertainties are 
given. Where a signed value is quoted, this indicates the direction in which as{44 GeV) 
changed with respect to the standard analysis. The scale uncertainty and quark mass 
effects are treated as asymmetric uncertainties of as- 
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1 -T 


Mh 


Bx 


Bw 


D2 


«,(35 GeV) 


0.1510 


0.1445 


0.1448 


0.1326 


0.1448 


fit range 


0.08-0.3 


0.22-0.46 


0.080-0.27 


0.06-0.16 


0.020-0.200 


xVd.o.f. 


25.2/8 


32.8/5 


23.7/8 


23.0/4 


20.3/8 


Statistical error 


±0.0009 


±0.0009 


±0.0007 


±0.0009 


±0.0014 


"I'Tcir'U'G r^nl^T^ 
licLL/iYO uiiiy 


— n nni (S 


<^ n nnm 


— n nm n 

U .UU lU 


— n nnnfi 

u.uuuu 


— n nm q 

U.UUlc/ 


LvlU-olt-io uiiiy 


-i-D nni 9 


-I-n nm ^ 


— n nnnn 

U.UUUc/ 


— n nm r 

U.UUIO 


— n nnnfi 

u.uuuu 




+n nnnzi 


-I-n nnn^ 


-I-n nnm 

znu.uuui 


-I-n nnn9 

znu.uuuz, 


-I-n nm 9 

znu.uui^ 


Pmiss 


-i-n nnm 


-I-n nnm 

nzu.uuui 


-I-n nnm 

nzu.uuui 


-I-n nnn'^ 

zizu.uuuo 


-I-n nnn4 

znu.uuui 


T)u 1 
ybal 


±0 0006 


±0 0001 


±0 0002 


±0 0006 


±0 0004 




+0.0006 


+0.0005 


+0.0005 


+0.0006 


+0.0005 




±0.0001 


±0.0001 


±0.0001 


±0.0001 


±0.0002 


fit range 


±0.0009 


±0.0017 


±0.0008 


±0.0016 


±0.0017 


Experimental syst. 


±0.0021 


±0.0024 


±0.0014 


±0.0026 


±0.0030 




-i-n nn98 


-I-n nn'?^ 

n^U.UUOtJ 


-I-n nn9'? 

n^U.UUZiO 


-I-n nm s 

n^u.uuio 


— n nnn9 

U.UUUZi 
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— u.uuoo 


— n nn9i 

— U.UUZl 


n nn9n 

— u.uu.^u 
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n^u.uuuo 


-I-n nm 

n^U .UU lO 
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n^U.UUltJ 
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n^u.uuuo 


rr -1- '^n A/TpV 
q I OU iVlc V 


— n nni ^ 


— n nnnfi 

u.uuuu 


— n nm 9 

U .UU IZi 


— n nm 

U.UUIO 


— n nnnzL 

u.uuuu 


1 .1 1 \ 1 I c^rin in p^'ti r* 
±jUiN-L' oV iliilieLi IL/ 


-I-n nnzin 


-I-n nn'^zL 

n^u.uuo^ 


-I-n nn97 

n^U.UUZi I 


-I-n nn9Q 

n^U.UUZiC/ 


-I-n nnnn 

n^U.UUUc/ 
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n nnnfi 


-I-n nm ^ 

n^U.UU lO 


— n nni4 

— U.UUl^ 
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T^U.UUl^ 


-un nn9zL 

T^U.UU.^^ 


ii^Q — ouu ivie V 


-un nnm 


n nnnfi 

— u .uuuu 


-un nnnfi 

n^u.uuuu 


— n nnns 

— U.UUUO 


— n nnnzL 

— u.uuuu 


A _ 50 MeV 


— n nnns 


— n noo6 


—0 0021 


— n noo3 


-)-n nnn9 


A + 50 MeV 


+0.0011 


+0.0007 


+0.0018 


+0.0003 


-0.0008 


udsc only 


+0.0074 


+0.0025 


+0.0086 


+0.0077 


+0.0055 


ML statistics 


±0.0008 


±0.0008 


±0.0007 


±0.0008 


±0.0013 


MC modelling 


+u.uuyz 
-0.0054 


+U.UUOU 
-0.0055 


+u.uuyy 
-0.0048 


+U.UUoy 
-0.0045 


+U.UUOO 
-0.0034 


= 0.5 


—0.0100 


—0.0077 


—0.0107 


—0.0078 


—0.0008 


= 2.0 


+0.0129 


+0.0103 


+0.0134 


+0.0097 


+0.0055 


Total error 


+0.0160 
-0.0116 


+0.0122 
-0.0098 


+0.0167 
-0.0118 


+0.0134 
-0.0094 


+0.0091 
-0.0048 



Table 8: Values of as(35 GeV) derived as in Table |^ but at 35 GeV. 
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D2 


a, (22 GeV) 


0.1607 


fit range 


0.060-0.200 


xVd.o.f. 


1.7/4 


Statistical error 


±0.0083 


tracks only 


+0 0023 


clusters only 


-0 0030 

KJ * KJ KJ KJ 


r iVTlPTirnprif J^l 
J_j.i\.lj\l;l lllidlljcxi oyoL. 


+0 0030 

_i_w . w wow 


a - 0.225 


-0.0015 


a + 0.225 


+0.0010 


ag - 30 MeV 


+0.0010 


(Tg + 30 MeV 


-0.0004 


LUND symmetric 


+0.0034 


Qo + 500 MeV 


+0.0031 


Qo - 500 MeV 


+0.0001 


A - 50 MeV 


+0.0011 


A + 50 MeV 


-0.0009 


udsc only 


+0.0105 


MC statistics 


±0.0025 


MC modelling 


+0.0119 
-0.0056 


= 0.5 


< 0.0001 


x^ = 2.0 


+0.0066 


Total error 


+0.0162 
-0.0105 



Table 9: Value of a;s(22 GeV) derived as in Table |^ but only for the differential 2-jet rate 
D2 at 22 GeV. 
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44 GeV 


1 - T 


Mu 


JJT' 


j-^ yy 


Do 


averaged 


a, (44 GeV) 


0.1510 


0.1532 


0.1681 


0.1406 


0.1302 


0.1442 


fit range 


0.12-0.35 


0.26-0.50 


0.16-0.35 


0.10-0.20 


0.01-0.20 




xVd.o.f. {x, = 1) 


3.0 


2.3 


3.7 


2.2 


1.1 




x^ ntted 


O.Ooo 


0.i3z 


0.600 


0.0 i\) 


O.OoO 




X^/d.o.f. (x^ free) 


2.0 


2.0 


3.6 


2.0 


1.8 




Statistical error 


±0.0028 


±0.0027 


±0.0025 


±0.0026 


±0.0021 


±0.0025 


Experimental syst. 


±0.0038 


±0.0051 


±0.0036 


±0.0034 


±0.0038 


±0.0029 


MC modelling 


±0.0041 
-0.0027 


±0.0034 
-0.0033 


±0.0074 
-0.0032 


±0.0040 
-0.0031 


±0.0057 
-0.0036 


±0.0048 
-0.0029 


Higher orders 


±0.0241 


±0.0117 


±0.0072 


±0.0086 


±0.0036 


±0.0072 


Total error 


±0.0249 
-0.0247 


±0.0135 
-0.0135 


±0.0112 
-0.0090 


±0.0104 
-0.0101 


±0.0080 
-0.0067 


±0.0095 
-0.0086 



Table 10: Values of as(44 GeV) derived using the 0{a1) QCD calculations with fixed 
x^ = 1 and fitted. The statistical and systematic uncertainties are also given. 



35 GeV 


1 -T 


Mh 


Bj' 


B\v 


D2 


averaged 


a, (35 GeV) 


0.1560 


0.1654 


0.1699 


0.1508 


0.1485 


0.1560 


fit range 


0.12-0.30 


0.30-0.50 


0.16-0.30 


0.10-0.20 


0.04-0.20 




xVd.o.f. {x, = 1) 


5.8 


2.3 


1.6 


3.4 


2.9 




fitted 


0.040 


0.342 


0.367 


0.056 


0.074 




X^/d.o.f. (x^ free) 


2.4 


2.2 


0.2 


1.4 


3.1 




Statistical error 


±0.0015 


±0.0019 


±0.0013 


±0.0014 


±0.0020 


±0.0016 


Experimental syst. 


±0.0034 


±0.0034 


±0.0025 


±0.0030 


±0.0052 


±0.0026 


MC modelling 


±0.0065 
-0.0036 


±0.0048 
-0.0048 


±0.0134 
-0.0064 


±0.0079 
-0.0047 


±0.0033 
-0.0032 


±0.0045 
-0.0029 


Higher orders 


±0.0279 


±0.0089 


±0.0140 


±0.0085 


±0.0068 


±0.0093 


Total error 


±0.0289 
-0.0284 


±0.0109 
-0.0108 


±0.0196 
-0.0157 


±0.0121 
-0.0103 


±0.0094 
-0.0094 


±0.0108 
-0.0102 



Table 11: Values of as(35 GeV) derived as in Table |lO|but at 35 GeV. 
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(a) 


(l-T) 


(Mys) 


(Br) 


(Bw) 


(y23) 


average 


aJMyo) 


0.1204 


0.1118 


0.1158 


0.1105 


0.1232 


0.1155 


Q range [GeV] 


13-172 


14-172 


35-161 


35-161 


22-161 




xVd.o.f. 


42.6/24 


10.9/14 


37.4/9 


21.1/9 


5.0/6 




experimental 


±0.0013 


±0.0010 


±0.0018 


±0.0015 


±0.0020 


±0.0013 


= 0.5 


-0.0050 


-0.0027 


-0.0039 


-0.0012 


-0.0043 


-0.0033 


x^ = 2.0 


+0.0061 


+0.0037 


+0.0048 


+0.0020 


+0.0057 


+0.0045 


fii = 1 GeV 


+0.0023 


+0.0021 


+0.0056 


+0.0054 




+0.0029 


Hi = 3 GeV 


-0.0018 


-0.0016 


-0.0043 


-0.0037 




-0.0019 


± 50% 




-0.0022 
+0.0025 


-0.0047 
+0.0064 


-0.0052 
+0.0073 




-0.0020 
+0.0028 


Total error 


+0.0066 

-0.0055 


+0.0050 
-0.0010 


+0.0099 
-0.0077 


+0.0094 
-0.0067 


+0.0060 
-0.0016 


+0.0062 
-0.0015 



(b) 


(l-T) 


(M-fj/s) 


(Br) 


(Bw) 




0.543 


0.457 


0.342 


0.264 


experimental 


±0.014 


±0.009 


±0.007 


±0.002 


X/^ — 0.5 


+0.002 


+0.013 


+0.009 


+0.030 


x^ = 2.0 


-0.001 


-0.008 


-0.006 


-0.019 


aj. ± 50% 




-0.076 
+0.212 


-0.036 
+0.063 


-0.024 
+0.037 


Total error 


+0.015 
-0.014 


+0.212 
-0.077 


+0.064 
-0.038 


+0.048 
-0.031 



Table 12: Values of Q;s(Mzo)(a) and «o (b) derived using the O^a"^) calculations and 
power corrections with iij = 2 GeV and x^ = 1. Fit ranges and x^/d.o.f. values for each of 
the five event shape observables are included. In addition, the statistical and systematic 
uncertainties are given. Where a signed value is quoted, this indicates the direction in 
which as{Mzo) and ao changed with respect to the standard analysis. The renormalisation 
and infrared scale uncertainties and the uncertainties due to the Ojf coefficients are treated 
as an asymmetric uncertainty on ^^(Mzo). These uncertainties are treated similarly for 
ao but exclude the infrared scale uncertainty. 
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Figures 





Figure 1: Measured and uncorrected distributions of the thrust observable 1— T (top) and 
of the differential 2-jet rate D2 (bottom) at 35 (left) and 44 GeV (right). The simulated 
data are overlayed as a solid line histogram. Only statistical errors are shown by the error 
bars. 
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Figure 2: The detector correction factors at 35 (points) and at 44 GeV (open squares) 
are sliown for the differential 2-jet rate, D2, in the Durham jet finder scheme. The error 
bars represent the statistical error. The arrow indicates the range of data considered to 
determine ag at 22 GeV. 
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Figure 3: Event shape distributions at \/s = 44 GeV corrected to the hadron level are 
shown for Thrust (T), heavy jet mass (Mh), total (-Bt) and wide jet broadening (Bw)- 
The error bars show the statistical error (inner tick marks) and the total error obtained 
by adding the statistical and experimental systematic error in quadrature. The solid line 
represents the JETSET 6.3 parton shower model prediction. 
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Figure 4: Event shape distributions corrected to the hadron level as for Figure || but at 
^/s = 35 GeV. The solid line represents the JETSET 6.3 parton shower model prediction. 
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Figure 5: Event shape distributions corrected to the hadron level at y/s = 44, 35, and 
22 GeV are shown for the differential 2-jet rate (-D2) in the Durham scheme. The error 
bars show the statistical error (inner tick marks) and the total error obtained by adding 
the statistical and experimental systematic error in quadrature. The solid line represents 
the JETSET 6.3 parton shower model prediction. 
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Figure 6: The distributions measured at ^/s = 44 GeV and corrected to parton level 
are shown for thrust T, heavy jet mass Mh, total and wide jet broadening Bt and Bw- 
The fits of the 0(a2)+NLLA (sohd line) and of the 0(af)(a;^ = 1) (dashed line) QCD 
predictions are overlayed and the fit ranges are indicated by the solid and dashed arrows. 
The error bars represent statistical errors only. 
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Figure 7: The same distributions as in Figure |] but for -^i = 35 GeV. 
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Figure 8: The distributions of the differential 2-jet rate, -D2, measured at ^/s = 44, 35, 
and 22 GeV using the Durham scheme are shown after correction to the parton leveL The 
sohd and dashed hues correspond to the fit results as in Figure 
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Figure 9: Values of as{35 GeV) and as(44 GeV) derived from 0(a^)+NLLA fits to event 
shape distributions. The experimental and statistical uncertainties are represented by the 
solid error bars. The dashed error bars show the total error including hadronisation and 
higher order effects. The shaded region shows the one standard deviation region around 
the weighted average (see text). For comparison the as values and errors measured by the 
OPAL Collaboration iSTI for the same set of observables are also shown. 
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Figure 10: Energy dependence of the mean values of thrust (1 — T), heavy jet mass 
{M]j/s), total (-Bt) and wide jet broadening {B\y), and of the differential 2-jet rate (7/23) 
are shown The solid curve is the result of the fit using perturbative calcu- 

lations plus power corrections while the dashed line is the perturbative prediction for the 
same value of ^^(Mzo). 
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Figure 11: Values of from (9(Q;f)+NLLA fits, as a function of centre-of-mass energy. 
The solid error bars are the statistical and experimental uncertainties added in quadrature, 
the dotted error bars are the total errors. The results from OPAL pn| , pT| for the same set 
of observables are shown as representative for the LEP experiments because the relevant 
detector subsystems of OPAL are similar to those of JADE. The solid line and the shaded 
band represent the QCD prediction for as{Mzo] 



0.1207±0.0012 corresponding to 



(242±15) MeV which was obtained from a x fit to the data taking only experimental errors 



into account 

shown by the dashed curves 



(5) 

For comparison the QCD predictions for A^ 



100 MeV and 400 MeV are 
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